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M(rjb) and ct s from Nonrelativistic Renormalization Group 

Bernd A. Kniehl, 1 Alexander A. Penin, 2,3 Antonio Pineda, 4 Vladimir A. Smirnov, 5 ' 1 Matthias Steinhauser 1 
1 77. Institut fur Theoretische Physik, Universitat Hamburg, 22761 Hamburg, Germany 
2 Institut fur Theoretische Teilchenphysik, Universitat Karlsruhe, 76128 Karlsruhe, Germany 

3 Institute for Nuclear Research, Russian Academy of Sciences, 117312 Moscow, Russia 
4 Dept. d'Estructura i Constituents de la Materia, U. Barcelona, E-08028 Barcelona, Spain 
5 Institute for Nuclear Physics, Moscow State University, 119899 Moscow, Russia 

(Dated: October 15, 2009) 

We sum up the next-to-leading logarithmic corrections to the heavy-quarkonium hyperfine split- 
ting using the nonrelativistic renormalization group. On the basis of this result, we predict the mass 
of the rjb meson to be M(r]b) = 9419 ± 11 (th) lg (5a s ) MeV. The experimental measurement of 
M(rjb) with a few MeV error would be sufficient to determine a s (Mz) with an accuracy of ±0.003. 
The use of the nonrelativistic renormalization group is mandatory to reproduce the experimental 
value of the hyperfine splitting in charmonium. 
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The theoretical study of nonrelativistic heavy-quark- 
antiquark systems is among the earliest applications of 
perturbative quantum chromodynamics (QCD) [l[ and 
has by now become a classical problem. Its appli- 
cations to bottomonium or toponium physics entirely 
rely on the first principles of QCD. This makes heavy- 
quark-antiquark systems an ideal laboratory to deter- 
mine fundamental parameters of QCD, such as the 
strong-coupling constant a s and the heavy-quark masses 
m q . Besides its phenomenological importance, the heavy- 
quarkonium system is also very interesting from the the- 
oretical point of view because it possesses a highly so- 
phisticated multiscale dynamics and its study demands 
the full power of the effective-field-theory approach. The 
properties of the T mesons, the bottom quark-antiquark 
spin-one bound states, are measured experimentally with 
great precision, and recent theoretical analysis of the T 
family based on high-order perturbative calculations re- 
sulted in determinations of the bottom-quark mass mt> 
with unprecedent accuracy @, H, 0] • 

In contrast to the T family, the current experimen- 
tal situation with the spin-zero rjb meson is rather un- 
certain: only one candidate event in 77 — > rjb produc- 
tion has been detected so far, which, however, is consis- 
tent with the expected background [5j. Yet, the discov- 
ery of the rj b meson is one of the primary goals of the 
CLEO-c research program Q. An accurate prediction 
of its mass M{jjb) is thus a big challenge and a test for 
the QCD theory of heavy quarkonium. Moreover, the 
hyperfine splitting (HFS) of the bottomonium ground 
state, E Ms = M(T(1S)) — M(r] b ), is very sensitive to 
a s and, with the advancement of the experimental mea- 
surements, could become a competitive source for the 
determination of the strong coupling constant. 

The HFS in quarkonium has been a subject of several 
theoretical researches Q. To our knowledge, the next-to- 
lcading order (NLO) G{a s ) correction is currently known 



in a closed analytical form only for the ground state HFS 
In this letter, we generalize this result to the excited 
states and present the analytical renormalization-group- 
improved expression for the heavy-quarkonium HFS in 
the next-to-leading logarithmic (NLL) approximation, 
which sums up all the corrections of the form a™ ln™^ 1 a s . 
We apply it to predict M{jjb). The result can be used for 
extracting a s (Mz) from future experimental data on the 
rjb meson mass. 

The leading-order (LO) result for the HFS is 
proportional to the fourth power of a s , E^® — 
Cf,a 4 s (ji)m q /{3n 3 ), where C F = (V c 2 - 1)/(2V C ), and 
suffers from a strong dependence on the renormalization 
scale /i of a s (/i), which essentially limits the numerical 
accuracy of the approximation. Thus, the proper fixing 
of ji is mandatory for the HFS phenomenology. The scale 
dependence of a finite-order result is canceled against 
the higher-order logarithmic contributions proportional 
to a power of ln(/i//i), where p, is a dynamical scale of 
the nonrelativistic bound-state problem. The physical 
choice of the scale fi = p, eliminates these potentially 
large logarithmic terms and a priori minimizes the scale 
dependence. However, the dynamics of the nonrelativis- 
tic bound state is characterized by three well separated 
scales: the hard scale of the heavy-quark mass m q , the 
soft scale of the bound-state momentum vm q , and the 
ultrasoft scale of the bound-state energy v 2 m q , where 
v oc a s is the velocity of the heavy quark inside the ap- 
proximately Coulombic bound state. To make the pro- 
cedure of scale fixing self-consistent, one has to resum to 
all orders the large logarithms of the scale ratios charac- 
teristic for the nonrelativistic bound-state problem. The 
resummation of the logarithmic corrections requires an 
appropriate conceptual framework. The effective field 
theory is now recognized as a powerful tool for the 
analysis of multiscale systems, which is at the heart of 
the recent progress in the perturbative QCD bound-state 
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calculations. The main idea of this method is to decom- 
pose the complicated multiscale problem into a sequence 
of simpler problems, each involving a smaller number of 
scales. The logarithmic corrections originate from loga- 
rithmic integrals over virtual momenta ranging between 
the scales and reveal themselves as the singularities of 
the effective-theory couplings. The renormalization of 
these singularities allows one to derive the equations of 
the nonrelativistic renormalization group (NRG), which 
describe the running of the effective-theory couplings, i.e. 
their dependence on the effective-theory cutoffs. The so- 
lution of these equations sums up the logarithms of the 
scale ratios. 

To derive the NRG equations necessary for the NLL 
analysis of the HFS, we rely on the method based on the 
formulation of the nonrelativistic effective theory known 
as potential NRQCD (pNRQCD) % The method was 
developed in Ref. [l(| where, in particular, the leading 
logarithmic (LL) result for the HFS has been obtained 
(see also Ref. [Hi). A characteristic feature of the NRG 
is the correlation of the dynamical scales, which leads to 
the correlation of the cutoffs jl2| ■ For perturbative calcu- 
lations within the effective theory, dimensional regular- 
ization is used to handle the divergences, and the formal 
expressions derived from the Feynman rules of the effec- 
tive theory are understood in the sense of the threshold 
expansion This approach 14, 15, H, 17 1 possesses 
two crucial virtues: the absence of additional regulator 
scales and the automatic matching of the contributions 
from different scales. 

Let us give a few details of the NLL analysis. We 
distinguish the soft, potential, and ultrasoft anomalous 
dimensions corresponding to the ultraviolet divergences 
of the soft, potential, and ultrasoft regions The LL 
approximation is determined by the one-loop soft run- 
ning of the effective Fermi coupling cf and the spin- flip 
four-quark operator [Ig| ■ In the NLL approximation, all 
three types of running contribute. We need the two- 
loop soft running of cp, which is known [l8j], and the 
two-loop soft running of the spin-flip four-quark opera- 
tor, which we compute by adopting the technique used 
in Ref. [3] for the calculation of the two-loop l/(m q r 2 ) 
non-Abelian potential. To compute the potential run- 
ning, we inspect all operators that lead to spin-dependent 
ultraviolet divergences in the time-independent pertur- 
bation th eory contribution with one and two potential 
loops [l(J[li|. They include (i) the 0(v 2 ,a s v) operators 
@, El, ( n ) the tree C(u 4 ) operators, some of which can 
be checked against the QED analysis [l6j], and (iii) the 
one-loop 0(a s v 3 ) operators, for which only the Abelian 
parts are known |l6[, while the non-Abelian parts are 
new. In the NLL approximation, we need the LL soft 
and ultrasoft running of the 0(v 2 ) and 0(v 4 ) operators, 
which enter the two-loop time-independent perturbation 
theory diagrams, and the NLL soft and ultrasoft running 
of the 0(a s v) and 0(a s v 3 ) operators, which contribute 



at one loop. The running of the Ojv 2 ,a s v) operators 
is already known within pNRQCD [10(. The running of 
the other operators is new. For some of them, it can be 
obtained using reparameterization invariance [2fJ. 

Besides the running discussed above, we need the ini- 
tial conditions for the NRG evolution given by the known 
one-loop result Q. With the anomalous dimensions and 
initial conditions at hand, it is straightforward to solve 
the system of the nonlinear differential equations for the 
effective couplings and get the NLL result for the HFS. 
The corresponding expression for general color (light fla- 
vor) number iV c (ni) and for arbitrary principal quantum 
number n is too lengthy to be shown in this Letter, so we 
present the explicit analytical expression only for N c = 3, 
ri; = 4, and n = 1, which applies to the bottomonium 
ground state. It reads 
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where a s is renormalized in the MS scheme, y = 
(Xsip) I OL a {rnb) , fl = CFOt s {[i)mb, B z (a,b) is the incom- 
plete beta-function, 2F\{a, b; c; z) is the hypergeometric 
function, and 2 -Fi(l, 1; 82/25; -1) = 0.7875078.... By 
expanding the resummed expression up to C(a 2 ), we get 
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where a s = a s (ju), fx = /x/n, CU = V c , T F = 1/2, 
/?i is the (i + l)-loop coefficient of the QCD (3 func- 
tion (/?o = HCa/3 - 4I>nj/3, . . .), L£ s = In (C F a a /n), 
*„(ar) = d"lnr(x)/dx™, r(x) is Euler's T function, and 
7s = 0.577216 ... is Euler's constant. In Eq. we keep 
the full dependence on N c , m, and n. The 0(a 2 In 2 a s ) 
term is known [I3,[lH, while the 0(a 2 In a s ) term is new. 

For the numerical estimates, we adopt the following 
strategy. First, we compute the dimensionless ratio 
Eh{ s /M(T(lS)), which only has logarithmic dependence 
on mb- At the order of interest, we have M(T(15)) = 
2m;,. Then, the ratio is converted into -E^fs by multiply- 
ing it with the physical T(15) meson mass, M (T(15)) = 
9460.30 ± 0.25 MeV 0. In Fig.HJ the HFS for the bot- 
tomonium ground state is plotted as a function of /i in the 
LO, NLO, LL, and NLL approximations. As we see, the 
LL curve shows a weaker scale dependence compared to 
the LO one. The scale dependence of the NLO and NLL 
expressions is further reduced, and, moreover, the NLL 
approximation remains stable at the physically motivated 
scale of the inverse Bohr radius, C F a s nib/2 ~ 1.5 GeV, 
where the fixed-order expansion breaks down. At the 
scale // w 1.3 GeV, which is close to the inverse Bohr 
radius, the NLL correction vanishes as one expects from 
general arguments. Furthermore, at //' w 1.4 GeV, the 
result becomes independent of fi; i.e., the NLL curve 
shows a local maximum. This suggests a nice conver- 
gence of the logarithmic expansion despite the presence 
of the ultrasoft contribution with a s normalized at the 
rather low scale p, 2 /nib ~ 0.8 GeV. By taking the differ- 
ence of the NLL and LL results at the local maxima as 
a conservative estimate of the error due to uncalculated 
higher-order contributions, we get -Ehfs = 41 ± 8 MeV. 

So far, we restricted the analysis to purely perturba- 
tive calculations. The nonperturbative contribution to 
the HFS can in principle be investigated by the method 
of vacuum condensate expansion 21]. The resulting se- 
ries does not converge well and suffers from large numer- 



ical uncertainties 22j. A reliable quantitative estimate 



of the nonperturbative contributions to the HFS can be 
obtained by comparison with lattice simulations. The 
HFS in bottomonium was studied on the lattice by sev- 
eral groups. The SESAM Collaboration 23j reported the 
value E his = 33.4 ± 1.9 MeV, while the CP-PACS Col- 
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FIG. 1: HFS of IS bottomonium as a function of the renor- 
malization scale fi in the LO (dotted line), NLO (dashed 
line), LL (dot-dashed line), and NLL (solid line) approxima- 
tions. For the NLL result, the band reflects the errors due to 
a s (M z ) =0.118 ±0.003. 



laboration [24| obtained E Ms = 33.2 ± 1.0 MeV (only 
the statistical errors are quoted). Both groups used un- 
quenched nonrelativistic lattice QCD. The result from 
quenched relativistic QCD on an anisotropic lattice [25| 
is -Ems = 33.8 ± 3.1 MeV. The central values are in good 
agreement and undershoot our NLL result by approxi- 
mately 7 MeV, which we take as an estimate of the non- 
perturbative contribution. 

Another way to estimate the accuracy of our result 
is by considering the HFS in charmonium, where ex- 
perimental data are available. The result of our anal- 
ysis is given in Figj2] along with the experimental value 
117.7 ± 1.3 MeV 0. The local maximum of the NLL 
curve corresponds to i?hfs = 112 MeV. Taking into ac- 
count that, in charmonium, the convergence of the per- 
turbative and relativistic expansion is much worse and 
the nonperturbative effects are much less suppressed in 
comparison to the bottom-quark case, the good agree- 
ment between our perturbative result and the experiment 
is amazing and supports our error estimates. We should 
emphasize the crucial role of the resummation to bring 
the perturbative analysis in agreement with the experi- 
mental data. Note that the recent lattice estimates un- 
dershoot the experimental value by 20 — 30% [271 ] . 

We may also apply our formulae to the n = 2 excited 
states. For bottomonium, a recent study of unquenched 
nonrelativistic lattice QCD with 1 + 2 light flavors by 
the HPQCD and UKQCD collaborations [H predicts 
the ratio E hls {2S) / E Ms (lS) = 0.32 ± 0.07, which is in 
good agreement with the value 0.27 obtained using our 
NLL result. For charmonium, our perturbative estimate 
Ehf s (2S) / Ehf s (lS) = 0.41 also agrees wery well with the 
result 0.41 ± 0.03 of the recent experimental measure- 
ments [29}. Although one definitely cannot rely on the 
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FIG. 2: HFS of IS charmonium as a function of the renor- 
malization scale /j, in the LO (dotted line), NLO (dashed 
line), LL (dot-dashed line), and NLL (solid line) approxima- 
tions. For the NLL result, the band reflects the errors due 
to a s (M z ) = 0.118 ± 0.003. The horizontal band gives the 
experimental value 117.7 ± 1.3 MeV 



(even NRG-improved) perturbative analysis of the ex- 
cited charmonium states, the above agreement suggests 
that the nonperturbative effects are not dominant and 
well under control at least for the ground state. 

In principle, our result can also be applied to the anal- 
ysis of the difference between the resonance energies in 
e + e~ — > ti and 77 — > ft threshold production [3fJ, which 
is beyond the scope of this Letter. 

To conclude, we have computed the heavy-quarkonium 
HFS in the NLL approximation by summing up the sub- 
leading logarithms of a s to all orders in the perturbative 
expansion. The use of the NRG extends the range of 
fi where the perturbative result is stable to the physical 
scale of the inverse Bohr radius. We found the resum- 
mation of logarithms to be crucial to reach agreement 
between our perturbative estimate and the experimental 
data on the HFS in charmonium despite a priori unsup- 
pressed nonperturbative effects. Our results further indi- 
cate that the properties of the physical charmonium and 
bottomonium ground states are dictated by perturbative 
dynamics. As an application of the result to the bot- 
tomonium spectrum, we predict the mass of the as-yet 
undiscovered 775 meson to be 
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